Interaction between particles in a 1/A-statistics anyon gas. 



(N 
O 
O 
(N 

C 



i 

CO 

-4-* 



I 

C 

o 
o 



> 
00 
On 

O 
(N 
O 

03 



o 
o 



X 
S3 



B. A. Bernevig + , D. Giuliano^, D. I. Santiago* 
+ Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 
t Universita di Napoli "Federico II" and I.N. F.M., Unitd di Napoli, Monte S.Angelo - via Cintia, 1-80126 Napoli, 
* Department of Physics, Stanford University, Stanford, California 94305 

(February 1, 2008) 



We analyze the dynamics between 1 / A- fractional statistics particles (anyons) in an exact three-body 
solution of the Sutherland Hamiltonian. We show that anyons interact by means of a short-range 
attraction. The interaction dictates important features of the spectral function for a charge- 1 
particle. 
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Particles constrained to move in a plane can exhibit 
exchange statistics intermediate between Fermi statistics 
and Bose statistics Intermediate-statistics particles 
are called anyons. Anyons occur in nature as the bulk 
excitations of Fractional Quantum Hall (FQH) samples 
0. In a FQH sample at filling fraction v — l/m, the 
elementary charged excitation is an anyon whose total 
charge is e/m and has l/m statistics. A similar frac- 
tionalization of electric charge and statistics takes place 
within the edge excitations of such samples. Therefore, 
edges of a FQH-sample provide a physical realization of 
a one-dimensional anyon gas |3) . Low-energy, long wave- 
length Quantum Hall edge excitations are described by 
a chiral Luttinger liquid model ||. The Luttinger liq- 
uid formalism, however, cannot provide a description of 
processes with high momentum transfer between anyons 
(short-distance physics). Short-distance anyon physics is 
what we study in this paper, in an exact solution of a 
simple one-dimensional model with 1/A-statistics excita- 
tions (A being a positive integer number): the Suther- 
land model H . The Hamiltonian for TV particles lying at 
xi, . . . , xn on a circle of length L and constant density 
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where Zi = exp[^p-]. 

By generalizing the formalism developed in || , we con- 
struct and analyze the Schrodinger equation for any num- 
ber of anyons. As a result, we show that anyons inter- 
act via a short-range attraction. When anyons are far 
apart, the probability of finding one at a certain position 
is essentially independent of the positions of the others, 
showing that anyons are free when far enough away from 
each other. When any number M (< A) of anyons come 



together there is an increase of probability. The maxi- 
mum enhancement of probability occurs when A of them 
come together in the 1/A statistics anyon gas. The short 
distance probability enhancement is the progeny of the 
short distance attraction; anyon attraction largely en- 
hances the probability of configurations with anyons at 
the same place. Although the contribution to the total 
energy due to the interaction disappears in the thermo- 
dynamic limit, the enhancement does not. Instead, the 
peak value of the probability increases with TV, and di- 
verges in the thermodynamic limit. The enhancement 
is the matrix element for a charge-1 excitation to decay 
into A- anyons. In the thermodynamic limit, the charge- 
1 excitation completely loses its integrity; as soon as a 
charge-1 particle/hole is created in the system, it imme- 
diately breaks up into A anyons. In particular for the 
A = 3 case we display the relevant physics of anyon at- 
traction in Fig.[j] where we plot the square modulus of the 
three-anyon wavefunction as a function of the separation 
between two anyons at a time. 
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FIG. 1. Square of the three-anyon wavefunction versus the 
separation between anyons at x a and xp (z-axis) and at xp 
and x-y (y-axis). Inset: spectral function vs. uj at q = 0.9n, 
q — 0.5-7T, q = 0. 

Let us start the mathematical derivation of anyon in- 
teraction. The ground state of Hs is 
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with energy £ = (^) A 2 ^|(iV 2 - 1). 

Elementary excitations on top of such a ground state 
are the one-dimensional analogs of Laughlin's anyons for 
a v = 1/A FQH- fluid. We create a 1/A charge and statis- 
tics quasihole at z a by removing 1/A of an electron from 
the ground state. The excitation is a localized charge 
defect given by the wavefunction: 

N 

^ a (zi, ■ ■ -,zn) = life' - z a )-$ GS (zi, . . . ,z N ) (3) 

3=1 

<J/ Q is not an eigenstate of Hs . Energy eigenstates are 
provided by Fourier transforming ty a : 

/dz 
^-^-^^ a { Zl ,...,z N ) (4) 

(to = 0,...,N and § denotes the integral over a 
closed path surrounding the origin). The correspond- 
ing eigenvalues are given by E + E(q m ) = E + 
(tt) 2 [(N - to) (Am + 1)]. The momentum of * m is de- 
fined as q m = (jf) [N/2 — to]. In the thermodynamic 
limit to/TV is kept finite, and E(q m ) — A [it 2 — , where 
-7T < q m < 7T. 

Multi-anyon states are constructed by a procedure 
analogous to the one used for one-anyon states. The 
wavefunction for M anyons localized at z ai , . . . , z aM is 
given by: 
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(5) 

A convenient basis for the M-anyon eigenstates is given 
by Fourier transforming Equation (|jj): 
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with < tom < . . . < mi < N. These many-anyon plane 
waves are not energy eigenstates. However, when acting 
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with Hs, we obtain 
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Eq.(^) is a lower-diagonal matrix equation. The en- 
ergy eigenvalues are the diagonal elements E{ m y — Eq + 

T,ili E (<l™,) - (w) Ei^i=i( m i - m i)- Many-anyon 
eigenstates of .ffj, <&{,„}, can be constructed from the 
^{m}' 13 by the so-called "squeezing procedure" Q. The 
(rrii—mj) terms make E^ m y not be equal to the sum of the 
energies of the single isolated anyons. Therefore, these 
terms show an interaction for anyons with large relative 
momenta, i.e at short distances. The anyon interaction 
is attractive because of the negative sign. 

In the thermodynamic limit, the energy eigenvalues 
E{ m y becomes equal to 
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We see that the energy for an M-anyon eigenstate just 
reduces to the sum of the energies of the single isolated 
anyons. The interaction contribution apparently disap- 
pears in the thermodynamic limit. However, as it hap- 
pens with spinons |j , the interaction plays a crucial role 
in determining some of the relevant physics of the anyon 
gas, as we are going to discuss next. 

The real-space wavefunction for M anyons localized 
at z ai , . . . , z au can be expanded in terms of the energy 
eigenstates: 
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where the function Vs m \ ( ^ 22 -, . . . , — — 
polynomial in each of its variables. By definition, 
nf=i C/^m} (fff, ■ ■ ■ , is the coordinate 

representation of the wavefuction for M anyons in the 
state of energy £{ m }. 

The Schrodinger equation for P{ m } is derived from the 
identity: 
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In Eq.(j7]) {to} is shorthand for 
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and from the fact that Hg acts on f ol! aM as: 
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a real electron out of the system. The spectral density 
of states, A(u, q) is the imaginary part of the Fourier 
transform of G with respect to both t and £, that is: 
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By putting together Eqs.(|9|, |10| , |ll[) , we obtain the equa- 
tion of motion for the relative coordinate M-anyon wave- 
function Vs m \ ( — , — , .... 
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(q = 2%n/L, < n < 3N is the total momentum of the 
hole). 

In order to compute A(uj,q), we need the norm of 
the energy eigenstates and the enhancements P{ m }(l, 1). 
The norms are obtained by following the methods devel- 
oped in H and can be expressed in terms of (^gsI^Gs) 
= T[3N + l]/T[3] N §. The solution is given by: 
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Eq.(|l2 

actly 



for three anyons and any A can be solved ex- 
|. In order to illustrate the features of anyon 
interaction, in the rest of this letter we shall analyze 
the full solution of Eq.([l2"|) for the case of three anyons 
and A = 3 (the case of two spinons, i.e. two anyons 
with A = 2 has been solved in ||.) [Pi m \(z, w)\ 2 , with 
z = z a2 /z ai w = z a3 /z a21 is the function we plot in 
Fig||, for A = 3, N = 20, m 1 -m 2 = 10, m 2 - m 3 = 7. 

We can extract physical consequences of anyon attrac- 
tion by computing the spectral density for the charge- 
1 excitation. The state where one more charge-1 hole 
has been created at £, on top of the TV-particle ground 
state, is described by the wavefunction ^(zi, . . . , zn) — 
n^=i( z i ~ 3i ^gs( z i, . . . , zn). From the definition of 
Vs m \ we find the following decomposition for ^ in terms 
of the states §{ m }'- 

*«= E E E r a+mp+m ^v {m} {i,i)^ {m} (is) 

m a —Q m^—Q m-y—0 

Since only positive-energy states propagate, the Green 
function for a charge-1 hole is defined as G(£, t) = 
-«0(*)(*£(f)|*£=i(t = 0)). The Green function corre- 
ponds to a physical hole, which is obtained by pulling 
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+ 1 + (j - i - l)/3]r[mj - m + ( j - i + l)/3] 



^ , {m}(ljl) can be computed from equation (|12[). The 
derivation will be described in a forthcoming paper [[|. 
Here, let us just quote the main result: 
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By inserting in Eq.(^) the formulas for the norms and 
for the enhancements, it is straightforward to work out 
A(u>, q) in the thermodynamic limit. In this limit, we can 
trade the sums for integrals over the Brillouin zone of the 
three anyons, so that Eq.(|l4|) takes the form: 



A(uJ,q) 
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In the inset of Fig.0, we plot A(q, u>) at different values 
of q. We recognize a fundamental feature of the anyon gas 
in that no definite Landau quasiparticle peak appears. 
This corresponds to the lack of integrity of the charge-1 
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hole we discussed before. More importantly, we can dis- 
tinguish a well-defined sharp peak in the spectral func- 
tion, whose height crucially depends on q. Such a "spec- 
tral enhancement" is the enhancement in the matrix ele- 
ment for a hole to decay in three anyons, V{ m j (1, 1). The 
dependence on q reflects the dependence of Vr m \(l, 1) on 
the relative anyon momenta, the larger the relative mo- 
menta, the larger the enhancement. Therefore, the max- 
imum enhancement should appear for values of q that 
make the differences among anyon momenta as large as 
possible, namely, when two anyons lie at opposite sites of 
the Brillouin zone (q « ±7r), while a third one has mo- 
mentum q « 0. On the other hand, the double constraint 
of momentum and energy conservation implies that the 
number of three-anyon states in which the hole decays is 
proportional to the size of the intersection between the 
plane J2i Qi^Q = 0; t ne sphere J2i ?f = 37r 2 — ui/6 and the 
Brillouin zone [— w, ir} 3 . Such a size is maximum when all 
three anyons lie at the corners of the Brillouin zone, that 
is, qi ~ ±7r i = 1,2, 3. The two competing maximum con- 
ditions discussed above imply that the maximum in the 
spectral function in this system does not appear exactly 
at the emission threshold, as it happens with spinons || , 
but at some intermediate value of u, depending on q. 

Such a trend clearly appears from the plots in the in- 
set of Fig|| as a signature of the nonelementarity of the 
charge- 1 excitation and of the corresponding attraction 
among fractionally charged quasiparticlcs. 

In the present work, the results of jBj have been gener- 
alized to excitations of arbitrary statistics 1/A. The main 
result is that fractionalizcd particles of any statistics al- 
ways interact through a short distance attraction discov- 
ered in H . We conclude that short distance attraction is 
a generic feature of the dynamics of fractionalized excita- 
tions. Such an attraction has profound consequences in 
the spectral function of the charge- 1 hole. The peculiar 
distribution of spectral weight should provide a possible 
experimental evidence of anyon attraction. 

In conclusion, it clearly appears that, because of its 
generality, our result should be a common feature of any 
anyon gas and, in particular, that it might provide a sim- 
ple explanation of some of the interesting and puzzling 
tunelling experimental findings on FQH samples [^],^| . 

As a mathematical appendix we provide the explicit 
recursion relation for the coefficients of V{ m y We assume 
a polynomial solution of the form 
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(15) 

where, if either p or q (or both of them) is < 0, b Piq = 0, 
bfi+kj = for k > j, & M+ „+i,« = for i > 0. 
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where the summations over p and q stop at some finite 
indices. The boundary condition implies feo,o = 1- More- 
over, by inserting into Eq.([l2]) the polynomial expansion, 
it is straightforward to work out a recursion relation for 
the coefficients b v „: 
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